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O ' Abstract 

In [Q], we present a method to tackle first order ordinary differential equations 
whose solutions contain Liouvillian functions (LFOODEs), many of them missed by 
QQ ■ the usual PS-approach. Here, we demonstrate an important result concerning the 

^^ ! generality of that method. 

O 

o 
o 

-(— > 



X 



*Umversidade do Estado do Rio de Janeiro, Instituto de Fisica, Departamento de Fisica Teorica, 
R. Sao Francisco Xavier, 524, Maracana, CEP 20550-013, Rio de Janeiro, RJ, Brazil. E-mail: 
lduarte@dft.if.uerj.br 

^idem. E-mail: sduarte@dft.if.uerj.br 

■'■idem. E-mail: damota@dft.if.uerj.br 

1 



1 Introduction 

The problem of solving ordinary differential equations (ODEs) has led, over the 
years, to a wide range of different methods for their solution. Along with the 
many techniques for calculating tricky integrals, these often occupy a large part of 
the mathematics syllabuses of university courses in applied mathematics round the 
world. 

The overwhelming majority of these methods are based on classification of the 
DE into types for which a method of solution is known, which has resulted in a gamut 
of methods that deal with specific classes of DEs. This scene changed somewhat 
at the end of the 19th century when Sophus Lie developed a general method to 
solve (or at least reduce the order of) ordinary differential equations (ODEs) given 
their symmetry transformations [§, y, Q. Lie's method is very powerful and highly 
general, but first requires that we find the symmetries of the differential equation, 
which may not be easy to do. Search methods have been developed fs], y] to extract 
the symmetries of a given ODE, however these methods are heuristic and cannot 
guarantee that, if symmetries exist, they will be found. 

A big step forward in constructing an algorithm for solving first order ODEs 
(FOODEs) analytically was taken in a seminal paper by Prelle and Singer (PS) [^ 
on autonomous systems of ODEs. Prelle and Singer's problem is equivalent to asking 
when a FOODE of the form y' = M(x,y)/N(x,y), with M and A^ polynomials 
in their arguments, has an elementary solution (a solution which can be written 
in terms of a combination of polynomials, logarithms, exponentials and radicals). 
Prelle and Singer were not exactly able to construct an algorithm for solving their 
problem, since they were not able to define a degree bound for the polynomials 
which might enter into the solution. Though this is important from a theoretical 
point of view, any pratical use of the PS method will have a degree bound imposed 
by the time necessary to perform the actual calculation needed to handle the ever- 
more complex equations. With this in mind it is possible to say that Prelle and 
Singer's original method is almost an algorithm, awaiting a theoretical degree bound 
to turn it algorithmic. In [|l|, we present a method which is an extension to the 
Prelle-Singer procedure allowing for the solution of some LFOODEs. This method 
is very effective in solving a class of LFOODEs missed by the usual Prelle-Singer 
procedure. It deppended though, for being general, on a result which we did not 
prove at the time. Here, we present such a proof thus guaranteeing that the method 
is generally applicable. The paper is organized as follows: in section g, we present a 
short theoretical introduction to the PS approach and show why the Prelle-Singer 
aproach misses some LFOODEs; in section ^, we demonstrate a theorem concerning 
the general structure of the integrating factor for LFOODEs of a certain type; we 
finally present our conclusions. 

2 The PS-Method and FOODEs with Liouvil- 
lian Functions in the Solutions 

2.1 The usual PS-Method 

Despite its usefulness in solving FOODEs, the Prelle-Singer procedure is not very 
well known outside mathematical circles, and so we present a brief overview of the 



main ideas of the procedure. Consider the class of FOODEs which can be written 

as 

I ^dy_ ^ M{x,y) 

^ dx N{x,y) ^ ^ 

where M(x,y) and N{x,y) are polynomials with coefficients in the complex field 
C. If i? is a integrating factor of (Q), it is possible to write dx{RN) + dy{RM) = 0, 
leading to N dr,R + R dyN + M dyR + R dyM = 0. Thus, we finally obtain: 

^ = -{d^N + dyM), (2) 

where D = Ndx + Mdy. In ||^], Prelle and Singer showed that, if the solution of ([l|) 
is written in terms of elementary functions, then R must be of the form R = Y\i /"' 
where fi are irreducible polynomials and rii are non-zero rational numbers. Using 
this result in (^), we have 



R n./r a/r 

i Ji' i J"- 

From (|2|), plus the fact that M and N are polynomials, we conclude that D[R]/R 
is a polynomial. One can then prove from (|3|) that fi\D[fi] [^. We now have a 
criterion for choosing the possible fi (build all the possible divisors of D[fi]) and, 
by using (|2|) and (^), we have 

Y.'^^ = - {d,N + dyM) . (4) 

i ■'» 

If we manage to solve (Q) and thereby find nj, we know the integrating factor for 
the FOODE and the problem is reduced to a quadrature. 

2.2 FOODEs with Liouvillian Functions in the Solu- 
tions 

In this section, we will explain why the PS-method fails to solve some FOODEs 
which solutions present Liouvillian functions. As we have already mentioned, the 
usual PS-method garantees to find a solution for a FOODE if it is expressible in 
terms of elementary functions. However, the method can also solve some FOODEs 
with non-elementary solutions. Why and why not for all the FOODEs with non- 
elementary solutions? Consider the following two examples: 

dy _ 3x2y2 + 2:^ + 1 

dx ~ 4(x + l) (x2 -x + l)y ^^ 

and 

-^ = y'^ + yx + x-l, (6) 



equation 1.18 of the standard testing ground for ODE solvers by Kamke M. These 
equations present (respectively) general solutions given by: 



y2 - V^TTVx2 -x + 1 1/2 / __ dx - 2_C1 = 0, (7) 

\ J Vx + \\Jx^ — x + 1 ) 

and 

— (8) 



2 _Ci + ^/^^e-2y2er/(l/2 V^^/2x - ^/^\/2) 

We note that both solutions are not expressible in terms of elementary functions. 
But, for FOODE ^, the standard PS-method can find the solution (eq. |^). The 
same is not true for FOODE y, what is happening? This can be best understood if 
we have a look on the integrating factors for those FOODEs, which are respectively: 

i2= (x=' + l)~^^^ (9) 

and 

x^ l2-2x 

R= ^. (10) 

(y + if 

Since the standard PS procedure constructs integrating factor candidates from poly- 
nomials in the variables {x,y), one can see that, since the integrating factor on eq. 
IC presents the exponential e^ /2-2x^ j^ ^jj^ never be found by the PS-method. 



So, can we understand something about the general structure of the FOODE that 
will allow us, eventually, to solve a class of these latter equations? In [|l|, we have 
demonstrated that, if (|l|) is a LFOODE, its integrating factor will be of the form: 

n 

R = e'°^^^y^Y[p,ix,yr^. (11) 

where tq is a rational function on {x,y), the pj's are irreducible polynomials on 
(x,y) and the Cj's are constants. Applying this into eq. (|2|), we get: 

D[roix, y)] + E ^^^^ = - {d.N + dyM) , (12) 

With this result, we could build a procedure to tackle certain classes of LFOODEs 
(missed by the usual PS-method) [^. In the bulk of the procedure, we use the 
hypothesis (based on our experience in dealing with those FOODEs) that pi are 
"eigenpolynomials" of the D operator. Our procedure proved to be successful in 
solving a class of such FOODEs and we actualy could not find a single counterexam- 
ple where our hypothesis (and consequently, our method) is not applicable. In the 
next section, we are going to prove that such a hypothesis is true and, consequently, 
our procedure is, indeed, generaly applicable. 

3 A Theorem concerning the structure of the 
integrating factor for a LFOODE 

As we have already mentioned, in ||l| we demonstrated that, for a LFOODE, the in- 
tegrating factor has the general form given by eq. (pT|) . Also there, we have prsented 



a method of solving LFOODEs that, with the (by the time of writing ||T|) hypothesis 
concerning the pj's (mentioned above), can be broadly divided into two steps: 

1. The determination of the pi's: Since they are "eigenpolynomials" of the 
D operator, is done in the same fashion as it is done in the usual PS-method 



(see section 2.1 above 



2. The determination of ro{x,y): In step 1, we get the pi's and, from eq.(12), 
we solve for ro{x,y), see [^. 

In this section, we are going to demonstrate that the above mentioned hypothesis 
is actualy a reality, therefore completing our analysis of the general form for the 
integrating factor for LFOODEs and, furthermore, assuring (by this demonstration) 
that the method presented on PI is generally applicable. Theorem 1: If we have 



a LFOODE of the form dy/dx = M{x,y)/N{x,y), where M{x,y) and N(x,y) are 
polynomials in {x,y), with integrating factor R given hy R = e^^^^'^' niLiK(3;,y)'^S 
where tq is rational function on (x, y), pi are irreducible polynomials on (x, y) and q 
are constants, then -D[ro] is a polynomial on (x, y), where D = Ndx + Mdy. Proof: 

Since tq is a rational function, we can write (|l^) as: 
'P{x,y) 



D 



Q{x,y) 



+ Y.Ci^^ = -id,N + dyM), (13) 



where P and Q are polynomials on (x, y) with no common factors. Writing J2i Ci—^ 
as a single quotient, we get: 



1 li Pi 



D 

Expanding D q and multiplying both sides of ( |l4|) by HiPj) we can write: 



n^^ ^^^^Q2^^^^^ + E c.( n P'i)D[P3\ = - id.N + dyM) (Up,). (15) 



Since X]j (^j(Wi,i^jPi)D\Pi\ is a polynomial and so is the right-hand side of (|15|) , we 
can conclude that 

i ^ 

is a polynomial. Since pi are irreducible polynomials, JJiPi can not divide Q^. So, 
we have two main situations: 

1. YliPi and Q have no common factors. 

2. YliPi and Q have common factors. 

Let us consider the first situation: Since, by assumption, Yl^pi does not have any 
common factor with Q and that Hi Pi — '~W — ^^ ^ polynomial, we must have 

that 

QD[P]-PD[Q] 

Q2 ^'') 



is itself a polynomial. Since {Q D[P] — P D[Q])/Q'^ = D[ro], the proof is complete 
for this case. The situation for the second case is a little more involved. First, let 
us consider that, in Hi Pi) ^ runs from 1 to n. With that in mind, let us stablish 
some notation: Representing the common factor of Q and HiLiPi as: 



n.i 



I=]Ip^, (18) 



i=l 



and the terms in 117= iP« ^^^ present in Q as: 



n 



^= U Pi^ (19) 

i=ni+l 

we can write: 

n 

l[p, = 7rl, Q = ei (20) 

j=i 

where 9 = Hl^i qT^ and rrij are positive integers. Re-writing (|lq ) with this notation 
and expanding, we obtain: 

^ QD[P]-PD[Q] QD[P]-PD[Q] 
1/^^ Q^ ' ^"^) Q9I 

QD[P]-PD[Q] D[P] D[Q] 

^" QO ^ "^-"^^F ^'') 

Remembering that ( ^l[ ) is a polynomial (see (jl^)), if we multiply ( ^I| ) by 9 (itself a 
polynomial), we get that 

ttD[P]-ttP^^ (22) 

is a polynomial. Therefore, since n D[P] is a polynomial, we finaly may conclude 
that 

.P^ (23) 

is a polynomial. From the fact that neither vr nor P have factors in common with 
Q, we can assure that D[Q]/Q is a polynomial. Let us consider the situation 2 
(therefore coompleting the possible cases): Reminding the reader that, by (pO|), 
Q = 9 I = (Yllli q^') iJXiliPi) and, using a similar line of reasonig that we used 
for R arround eq.(^), we can conclude that qi\D[qi\,i = l..ng and pi\D\pi\,i = l..ni. 
Writing eq.(12) as: 



D[ro{x,y)]+f:'-^^^+ ± '^^^^ = - (d.N + dyM) , (24) 

i=l P^ i=ni+l P' 

and, since pi\D[pi],i = l..ni, we have that J27=i ^' . is a polynomial. Therefore 

ON+ t aajll (25) 

j=ni+l ^« 



is a polynomial. We can write ( ^51) as: 



riiK 



6 



Multiplying (EBI) by UiPi, we get: 



Y[p^Q3n^^M + Y^c,{]lp,)D[p,], ^,J=n, + l...n (27) 

which is also a polynomial. Since J2j '^j(Jli,ij^j Pi)D[Pj] is itself a polynomial, we 
can safely say that: 

j=ni+l ^ 

is a polynomial. Since OiLni+i Pi has no common factor with Q, we finally conclude 
that — — —q2 — — = D[ro] is polynomial. Thus, we have demonstrated that D[ro\ 
is always polynomial. Corollary 1: If R = e^°^^'y'YVi^iPi{x,yy^ (where ro is a 



rational function on (x, y), the pi's are irreducible polynomials on (x, y) and the Cj's 
are constants) is the integrating factor for the LFOODE y' = M/N, where M, A'" 
are polynomials in (x,y), then pi\D[pi]. The result above is a direct consequence of 
eq. (PI). 



4 Conclusion 

In |l|], we have developed a method to deal with a class of LFOODEs. This method 
proved to be very effective and we could not find (within the class under study) 
any counterexample where it was not applicable. Nevertheless, in order to prove its 
genetrality, the result presented on this paper was lacking. For that reason only, 
our present result is justified. But, apart from that, it is worth mentioning that our 
result completes the determination of the general form for the integrating factor 
for any LFOODE of the type given by eq. ([l|). Now we can safelly say that the 
integrating factor for any LFOODE has the form: 

i? = e'■«(^'^')^p,(x,y)=^ (29) 

i=l 

where -D[ro] is a rational function on {x,y), the pi's are "eigenpolynomials" of the 
D operator and the Cj's are constants. 
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